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Abstract 

A 5C/(4) flavor symmetrical Lagrangian is constructed for the interaction of the 
pseudo-scalar mesons with the vector mesons. SU{4:) symmetry is broken to SU{3) 
by suppression of terms in the Lagrangian where the interaction should be driven 
by charmed mesons. Chiral symmetry can be restored by setting this new SU{4) 
symmetry breaking parameters to zero. Unitarization in coupled channels leads to the 
dynamical generation of resonances. Many known axial resonances can be identified 
including the new controversial A (3872) and the structure found recently by Belle 
around 3875 MeV in the hidden charm sector. Also new resonances are predicted, 
some of them with exotic quantum numbers. 

1 Introduction 

Recent years have been very exciting for the hadron spectroscopy because of the 
discovery of many new and controversial states that do not fit well the interpretation 
of baryons as qqq states or mesons as qq states. In particular two charmed resonances 
discovered by BaBar [T] and confirmed by other experiments [2], [3], [4], the Dso(2317) 
and I?si(2460) have animated the debate about non qq mesons. Also non-strange 
partners of these resonances have been observed [5], [6]. 

The predictions for the masses of these states with quark model potentials al- 
ready existed [7] and turned out to be off by more than 100 MeV. The fact that 
the Dso(2317) lies just below the DK threshold and the -Dsi(2460) just below the 
D*K threshold made many theoreticians speculate that these states could be meson 
molecules [8], [9], [lO], [H], [l2], [13]. Others support a tetraquark assignment [H], 
|15| . |16| . or usual qq states with more sophisticated quark model potentials or within 
QCD sum rules calculations [17], [18], [H], |20j and there is also the possibility of 
admixture between these configurations [21], [22], [23]. 

In the hidden-charm sector also new controversial resonances have been found. In 
particular the A (3872), observed in four different experiments [23], [25], [26], [27], has 
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attracted much attention. The narrow width of this state makes its interpretation 
as a usual charmonium cc state very difficult. For this resonance too, many exotic 
theoretical interpretations have been investigated such as tetraquarks, hybrids and 
molecules [28] , [29] , [30] , [31] , [32] . For a good review on heavy mesons one can refer 
to [33]. 

In this work unitarization in coupled channels is used to explore the pseudo-scalar 
meson interaction with vector mesons. In the works of Kolomeitsev [9] and Guo [H] a 
similar approach has been done using a Lagrangian based on heavy quark chiral sym- 
metry that allowed the investigation only of the open-charm sector and constrained 
the interaction for only light pseudo-scalars with heavy vector mesons. For our phe- 
nomenological model we construct a Lagrangian based on 5f7(4) flavor symmetry 
and this symmetry is broken to SU{3) by suppressing exchanges of heavy mesons 
in the implicit Weinberg- Tomozawa term. Chiral symmetry can be restored from 
our model by setting the S'C/(4) symmetry breaking parameters to zero. This new 
Lagrangian, based on the ideas of a previous paper [34], includes also the possibility 
to investigate the hidden-charm sector and the interaction of heavy pseudo-scalars 
with light vector mesons, which enriches the spectrum of axial resonances generated. 

The paper is organized as follows: in the next section the construction of the 
Lagrangian is explained in detail and also the mathematical framework for solving 
the scattering equations in a unitarized approach is presented. Section 3 is devoted 
to the presentation and discussion of the results and section 4 contemplates overview 
and conclusions. 



2 Mathematical Framework 

We will start by constructing two fields, one for the 5'C/(4) 15-plet of pseudo-scalars 
and another one for the 15-plet of vector mesons: 
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Now for each one of these fields a vector current is build: 



= (9^$)$- $5^$ (3) 

= (^^.v,)^^ - v.a^v^ (4) 

The Lagrangian is then constructed by connecting the two currents: 

^=lj2TriJ,Jn- (5) 

Note that the lo appearing in eq. ([2]) is not the physical lj but tog. The addition 
of a singlet state, which is diagonal and proportional to the identity matrix in the 
representation of eq. ([2]) does not give any contribution to J^^ in eq. (jj]) and hence 
does not modify the Lagrangian in eq. ([5]). However, there is a caveat since the singlet 
and octet mix strongly to give the uj and (j) states which have different masses. In 
the results section we shall come back to this problem and will discuss the effect of 
this mixing. 

Next step is to break SU (4) symmetry, this will be done by suppressing exchanges 
of heavy mesons. To identify these terms, one should first decompose each one of the 
fields into its SU{3) components: 
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The I3 is the 3x3 identity matrix and the fields (pi and contain the meson fields 
for each i-plet of SU (3) into which the 15-plet of S'C/(4) decomposes: 
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In terms of the 5f7(3) fields the Lagrangian reads: 



C = 



2 4 \ 

-^{JsSf^J^l + Jl3fiJQ + -^Jl3^lJ^lj + JsSf.J'^s + '^38m^8'3 + 

2 4 \ 

"y^(-'38/iv7i3 + -^31/i^83) + 2'^3lAf^l3j- 

In eq.(l8|) the currents are defined as: 
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The interaction in the Lagrangian of eq. ([5]) is usually visualized, in the vector- 
meson dominance picture, as a t-channel exchange of a vector meson in between 
the meson pairs. Since this Lagrangian is 5C/(4) flavor symmetric, it assumes equal 
masses for the virtual vector meson exchanged in between the meson pairs. In order 
to break SC/(4) symmetry we will suppress the terms in this Lagrangian where the 
virtual vector meson exchanged is a heavy one. If the transferred momentum in 
between the meson pairs is neglected, the propagator of the exchanged meson is 
proportional to the inverse of its squared mass {.^t\ Considering the terms where a 
light vector mesons is exchanged of the order of magnitude of the unit, the suppressed 
terms should be thought to be of the order where and mii are scales of 

the order of magnitude of the light and heavy vector mesons masses respectively. 

Terms in the Lagrangian where the two currents of eqs. ([9]) and (jlOh have ex- 
plicit charm quantum number should be suppressed since these currents can only 
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be connected through the exchange of a charmed, and hence heavy, vector meson. 
There are stiU two terms in the Lagrangian where a hidden-charm meson can be 
exchanged in between the meson pairs. These are the terms involving only mesons 
belonging to the SU{3) triplet and anti-triplet. In this case the interaction is driven 
simultaneously by light and heavy vector mesons and in order to correctly suppress 
these terms one should isolate each one of these contributions and suppress only the 
one coming from the heavy vector meson {J/ip). These contributions have already 
been calculated in a previous paper |34j and they can be either | and | for light 
and heavy vector mesons respectively, if two equal currents are connected, or — | 
and I for light and heavy vector mesons respectively, if two different currents are 
connected. So, following the steps in [M], the corrected Lagrangian, accounting for 
the masses of the heavy vector mesons, reads as: 

■^(-'83^^31 + '^ISf^J'^s'l + Y'^ISm^s'i) + ^^-^33^^33 + 

7<^38Af^83 + ■^(•^38Af^l3 + JufiJ^^) + "^'^SlAt-^ls) ' i^^) 

-ith7=fe)'andV = -i + t(^)'. 

The parameters niL, mn and m'^ should be chosen of the order of magnitude of a 
light vector-meson mass, charmed vector-meson mass and the J ftp mass, respectively. 
The masses of the light vector mesons vary in between 770 MeV, for the p meson 
and 892 MeV for the K* mass, the charmed vector-mesons have masses 2008 MeV 
and 2112 MeV {D* and D*, respectively) and the J/ip mass is approximately 3097 
MeV. We will chose mL=800 MeV, mH=2050 MeV and 771'^=^ GeV. Changing these 
parameters over the whole allowed physical range (770 MeV < mi < 892 MeV and 
2008 MeV < mn < 2112 MeV) has about 0.1% effect over the pole position of the 
heavy resonances. 

Note that the first term of the Lagrangian in eq. (jlip is the same Lagrangian used 
in the works of Lutz [35] and Roca [36] in the study of the low lying axial mesons 
and the third term is the same Lagrangian used in the works of Kolomeitsev [9] and 
Guo [11] when studying the open-charm sector. This matching sets the value as 
the coefficient in the Lagrangian of eq. dSj). 

From the Lagrangian in eq. pT|) one gets the transition amplitudes between an 
initial and a final state: 



M'rjis,t,n) = ^(s-n)e.e'. (12) 



The super-index C refers to the charge basis, and the labels i and j to the initial 
and final channels while s, t and u are the usual Mandelstam variables. In appendix 
A we give tables for the coefficients in an isospin basis. 
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These amplitudes will first be projected in s-wave: 



V/j{s) = d{cose)Mij(^s,t{s,cose),u{s,cosd)^ (13) 

This potential will then be used as kernel in a Bethe-Salpeter equation, in an 
on-shell formalism [37] , [38] , [39] , [40] . In this way the unitary T-matrix assumes the 
form |36j : 

T = -{i + VdyWl'.-t' (14) 
In this equation G is a diagonal matrix with each element given by: 



^ 1 / mf Mf -mf + s M? 

= i6^r+^°v+ — — 

^ (log + + ^pV^ I Ioq ^ + ~ + )\ (IQ) 

VsV ^-s + M^-mf + 2p^ ^-s-M^ + mf + 2p^J J ^ ' 

in the above equations p is the three-momentum in the center of mass frame of the 
two mesons in channel I, while Mi and m; are the masses of the vector and pseudo- 
scalar mesons respectively, a is the subtraction constant, which should be fitted as 
a free parameter and /i is a cut-off scale which we will set to 1.5 GeV. Note that a 
and II are not independent, this justifies setting one to a fixed value and adjusting 
just the other one to data. 

The three-momentum p is calculated with: 



^ = ^ ^ 

Unitarity is ensured by the imaginary part of the loop function of equation (jl6p : 

Im{Gi) = --^ (18) 

When looking for poles in the complex plane one should be careful because of 
the cuts of the loop function beyond each threshold. Bound states appear as poles 
over the real axis and below threshold in the first Riemann sheet. Resonances show 
themselves as poles above threshold and in the second Riemann sheet of the channels 
which are open. 

Over the real axis the discontinuity of the loop function is known to be two 
times its imaginary part [JT] so, knowing the value of the imaginary part of the loop 
function over the axis, eq. (jlSp . one can do a proper analytic continuation of it for 
the whole complex plane: 
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Gj' = + lm{p)>0 (19) 

G^^ and G^ refer to the loop function in the second and first Riemann sheets, respec- 
tively. 

Until now our formalism worked only with stable particles, but in some cases one 
has a /? or a K* meson in the coupled channels, and these particles have relatively 
large widths. The consideration of the mass distributions of these particles can be 
relevant whenever thresholds are open thanks to this mass distribution. 

In order to take this into account we follow the procedure of [32] and convolute 
the loop function with the spectral function of the particle, hence, using a new loop 
function: 



~ _ 1 r{MR+2rR)^ ~o/-l\ 1 

G{^,m,MR) = - dM^ — )lm^ 

N J{M„-2T„)'i V vr y M 



l{MR-2rnV V vr y M2 - M2 + iMnVR 

G{^s,m,M) (20) 

N = / dM^ — ■ (21) 



(Mfi-2r«)2 V vr y M2 - M2 + iMuTR 

In the next section we will comment further on this issue and present results, for 
the heavy resonances, by taking into account the finite width of the p and K* vector 
mesons in the few cases where the generated resonances have important coupling to 
channels involving these mesons and a mass close to the threshold of these channels. 



3 Results 

The 15-plet of SU{A) breaks down into 4 multiplets of 5C/(3): 

15 — > 1030308. (22) 

Knowing this, one can study the SU{3) structure of the interaction between 
pseudo-scalar and vector mesons. Table [D shows the SU{3) decomposition of the 
interaction. The irreducible representations (irreps) marked with a * refer to the 
vector meson multiplet. 

Since now one can differentiate between the vector and the pseudo-scalar repre- 
sentations in the irrep products, the coupled channel space gets enlarged with respect 
to the case of the scalar resonances and therefore, a much richer spectrum is gen- 
erated for the axial resonances. In contrast with previous works of Kolomeitsev [9] 
and Guo [11] where light pseudo-scalar mesons are scattered off heavy vector mesons, 
in our work we have an enlarged coupled channel basis accounting also for channels 
where light vector mesons are scattered off heavy pseudo-scalars. 

Using SU{3) isoscalar factors [34], [44], [45] the ^{j can be transformed to a SU{3) 
basis, by means of which one knows in which multiplets there is attraction and, 
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Table 1: SU{3) decomposition of the interaction between pseudo-scalar and vector mesons 
in SU{A). The sectors not shown in the table correspond to the C = —1, —2 states which 
are just charge conjugate states (antiparticles) from the ones shown. 



therefore, the possibility to generate resonances. In the C=2 sector the interaction 
is not attractive in any multiplet, in the C=l sector there is attraction in the anti- 
triplets and in the sextets and in the C=0 sector the two octets and the singlet 
coming from 8 8 8* and the two heavy singlets from the 3 8 3* and the 3 (8) 3* are 
attractive. 

Table [2] shows the channel content in each sector. We will work with states of 
defined charge conjugation or G-parity, where it applies. 

Our model assumes that isospin symmetry is exact so, all particles belonging to 
a same isospin multiplet have the same mass. For the pions we use m7r=138 MeV, 
for kaons m/^=495 MeV and for the eta m^=548 MeV. In the heavy sector we use, 
for the pseudo-scalars mD=1865 MeV, m£)^=1968 MeV and mr,^=2979 MeV. The 
masses of the light vector mesons are: mp=771 MeV, mi^* =892 MeV and m^=7S2 
MeV. While for the heavy vector mesons we use: m£)*=2008 MeV, m£)*=2112 MeV 
and mj/^=3097 MeV. 

It is also possible to restore SU{3) symmetry by setting the masses of all particles 
in a same SU{3) multiplet to a common value. For this purpose we introduce the 
parameter x, x = is the case when SU{3) symmetry is restored and x = 1 the case 
we see in Nature with SU{3) broken. The meson masses as a function of x are given 
by: 

m{x) = fh + x{mphys. — ™) (23) 
where fh is the meson mass in the symmetric limit. We will use for the pseudo-scalars: 
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Table 2: Channel content in each sector 
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Figure 1: Pole trajectories for one of the anti-triplets in the C=l sector while breaking 
SU{?>) symmetry in steps of Ax=0.1. The two degenerate poles at x=0 become two 
different resonances at x=l, in the two extremes of the curve (real world). 



777,8=430 MeV, m3=1900 MeV and mi=1900 MeV, and for the vector mesons: 
m8*=800 MeV, ?fi3*=2050 MeV and mi.=2050 MeV 

In our model we still have to fit the subtraction constants in the loop function. As 
done in our previous work [3l] we use two different values, for channels involving 
just light mesons and an for channels involving at least one heavy meson. This is 
justified because the heavy and light sector nearly decouple from each other. The 
result of our fit is ul = —0.8 and an = —1.55. 

We also use two different values for the meson decay constant / appearing in the 
Lagrangian: for light mesons / = /tt = 93 MeV, while for heavy mesons f = fo = 165 
MeV [43]. 

When SU{2>) symmetry is restored it is possible to identify 3 poles and one cusp 
in the open charm sector. Two poles come from the two anti-triplets where the 
interaction was attractive, their positions are at 2432.63 MeV and (2535. 07-10.08) 
MeV. The other pole is broad and comes from one of the sextets at (2532. 57-1199.36) 
MeV, the other sextet appears as a narrow cusp around 2700 MeV, it becomes a pole 
when the heavy light threshold at 2700 MeV moves because of the SU{2>) symmetry 
breaking. 

Figure [1] shows the pole trajectories for the anti-triplet starting at 2432.63 MeV 
while changing x from to 1 in steps of 0.1. 

In the hidden charm sector two octets and three singlets are expected, one light 
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and two heavy. The two octets are nearly degenerated at 1161.06 MeV and 1161.37 
MeV. In the work of Roca [36] the two octets are degenerated, but in our model the 
interaction with the heavy sector removes this degeneracy by a very small amount, 
indicating a weak coupling between the light and heavy sectors. The light singlet 
appears as a pole at 1055.77 MeV and the two heavy ones at 3867.59 MeV and 
(3864.62-iO.OO) MeV this second one is not exactly a bound state as the others, but 
a narrow state with a width smaller than 1 KeV. 

Table [3] shows the pole positions for the case x = 1 within our model and the 
possible identification of each one. 

We will now discuss separately the particle identification in each sector. 



3.1 C=1,S=1,I=1 

In contrast with the scalar resonances where the sextet state became very broad [34], 
the axial sextets are narrower, hence easier to detect experimentally. One should 
note also that these states are truly exotics since quark models cannot generate qq 
pairs with such quantum numbers. We found two poles in this sector at positions 
(2529. 30-1238. 56) MeV and (2756.52-132. 95) MeV. 

The couplings, gi, of the poles to each channel i have been calculated from the 
residues of each pole. Close to the pole position one can write: 



T,^ ^ _Ji9i_ (24) 

Table m shows the results of gi for the poles in this sector. With the couplings it 
is possible to do a rough estimate of the partial decay widths for the resonances and 
thus identify the channels with largest contribution to the width in order to motivate 
experimental searches in this direction. At tree level one has: 



^^-^^ = 8^^ 

where p is the center of mass three-momentum of the two particles in the final state. 
For the pole at (2756. 52-132. 95) MeV the estimate reads: 

~ 3.3 and ~ 3.6 



While for the pole at (2529.30-1238.56) MeV on has: 

— ^ ~ 3.1 

^ KD* 

The large coupling of the lighter state to ttD* and KD* , or the heavier one to Dgp 
and DK* make these states qualify as roughly quasi-bound states of these channels 
respectively. Note that they separate two basic configurations: heavy vector-light 
pseudo-scalar and heavy pseudo-scalar-light vector. 

When taking into account the finite p and K* widths the resonance at (2756.52- 
132.95) becomes a little bit higher in mass, crossing the KD* thresholds and disap- 
pearing as a pole. 
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Mass (MeV) 
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Table 3: Pole positions for the model. The column Irrep shows the results in the SU{3) 
limit. The results in brackets for the Irriy/s are obtained taking into account the finite 
width of the p and K* mesons. 
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Channel 


(2529.30-i238.56) MeV 

h\ (GeV) 


(2756.52-i32.95) MeV 
h\ (GeV) 


ttDI 


9.13 


2.54 


DsP 


2.18 


9.26 


DK* 


1.94 


11.02 


KD* 


8.20 


2.61 



Table 4: Residues for the C=1,S= 1,1=1 sector 



Channel 


2455.91 MeV 


(2573.62-i0.07) MeV 




1^7.1 (GeV) 


1^7.1 (GeV) 


DK* 


0.54 


13.96 


KD* 


9.74 


0.30 




6.00 


0.18 


DsUJ 


0.51 


7.95 


vcd: 


0.02 


0.05 


DsJ/^ 


0.54 


0.00 



Table 5: Residues for the C=1,S=1,I=0 sector 



3.2 C=1,S=1,I=0 

The two poles found in this sector have the proper quantum numbers to be identified 
with the two Dsi resonances. The first pole appears as an exact bound state at 
2455 MeV and we identified it with the Dsi(2460) state. Experimentally the main 
hadronic decay channel for this resonance is D*Tr which is an isospin violating decay 
and therefore not taken into account by our model. Other decays for this resonance 
are three body decays or electromagnetic ones, which are also not included in our 
framework. 

The other pole appears at (2573. 62-10.07) MeV and couples mainly to the DK* 
and DgLV channels. The only open channel for it to decay is the KD* channel but, 
because of the dynamics of the interaction, this resonance barely couples to it. This 
explains the small width of this resonance, 140 KeV, despite the 70 MeV phase space 
available for it to decay. We identify this pole with the Dsi(2536) which is also 
observed in the decay channel KD* with a small width (T < 2.3 MeV |43|). 

Table [5] shows the absolute value of the couplings \gi\ for each channel for the two 
poles in this sector. 

Once more we see that the lighter state couples strongly to KD* and r]D* while the 
second one couples strongly to DK* and Dgto. Hence the decoupling into two families 
of heavy vector-light pseudo-scalar and light vector-heavy pseudo-scalar shows up in 
this sector too. 

This sector has one resonance with a strong coupling to the Dgco channel {uj = los 
here). Hence, this is one case where the (j) — oj mixing can be relevant. Since we saw 
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Channel 


2455.97 MeV 


(2597.44-i0.06) MeV 




1^7.1 (GeV) 


1^7.1 (GeV) 


DK* 


0.41 


13.86 


KD* 


9.82 


0.27 


vD*s 


6.06 


0.09 


DsUJ 


0.27 


5.12 


Ds(l) 


0.36 


7.06 


VcD: 


0.02 


0.01 




0.55 


0.00 



Table 6: Residues for the C=1,S=1,I=0 sector when introducing the Dgcf) channel 



that the singlet (wi) does not give any contribution to the Lagrangian, the exphcit 
introduction of the u and ^ states can be done by substituting: 

This splitting into two fields introduces a new column and a new row in the second 
table of appendix A.I., by simply multiplying the fourth column and row by and 

introducing a Ds4> column and row with weights those of the original DgUJ^. 

When we look now for poles we obtain the results in Table El 

Comparing Tables [5] and [6l we can see that in the case where the coupling of the 
resonance to the DgU^ channel is weak (first resonance) the effects of the mixing are 
very small in the energy and couplings of the resonance. In the case of the second 
resonance, where the coupling to the Dsio% channel was large, the effects of the mixing 
are more visible. There is a shift of the mass of about 25 MeV, which is well within 
our theoretical uncertainties. This effect can be considered an upper bound for all 
other cases, since we have chosen the resonance with strongest coupling to wg- 

It is also interesting to see that the sums of the squares of the couplings to to and 
(j) are close to the square of that to tijg, indicating a redistribution of the strength of 
the coupling to oj^ between uj and cj). 

The widths of the light vector mesons have no significant effects over the reso- 
nances generated in this sector, because the mass of the resonances are far away from 
the threshold of the DK* channel. 

3.3 C=l,S=0,I=i 

Here the companions of the two anti-triplets and the two sextets should be found. 
Note that when we refer to the SU{2>) multiplet we are talking about the case when 
one has S'f7(3) symmetry. This correspond to x = in the pole trajectories. At x = 1, 
since SU{3) symmetry is broken, the physical states mix the SU{3) multiplets. Yet, 
the study of the trajectories allows us to trace back any pole to its origin in the SU{3) 
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Channel 


(2311.24-ill5.68) MeV 


(2526.47-i0.08) MeV 


(2750.22-i99.91) MeV 




h\ (GeV) 


h\ (GeV) 


\g,\ (GeV) 


TlD* 


9.84 


0.24 


2.15 


Dp 


0.89 


12.13 


3.82 


KDl 


5.21 


0.59 


2.38 


DsK* 


0.09 


7.89 


13.11 


r]D* 


0.68 


0.56 


2.43 


Du 


0.61 


1.58 


7.47 


VcD* 


0.06 


0.02 


0.22 


DJ/iP 


1.27 


0.02 


0.01 



Table 7: Residues for the C=1,S=0,I=| sector 



symmetric case, and we have used this information for the classification of states in 
Table Ej 

The anti-triplet companion of the pole for the Z)si(2460) is the pole located at 
(2311. 24-1115.68) MeV that we identify with the L>i(2430) because of its naturally 
large width, since it is strongly coupled to the ttD* channel into which it is free to 
decay. On the other hand the pole at (2526. 47-10.08) MeV, companion of the one 
identified with the L>si(2536), has its coupling to the irD* channel strongly suppressed 
and therefore has a very narrow width. Because of this unnatural narrow width we 
are tempted to identify it with the Di (2420) although the mass of our dynamically 
generated resonance is around 100 MeV off the experimental value for this state. 
Moreover, when considering the finite widths of the p and K* mesons, this pole gets 
a larger width, its imaginary part goes to —13 MeV, implying a width of about 26 
MeV, in fair agreement with experiment. 

As for the sextets, one of the poles becomes a broad cusp at the KD* threshold 
as one gradually breaks SU{3) symmetry through the parameter x, and the other 
pole emerges from a cusp into a pole at (2750. 22-199. 91) MeV. The channel to which 
it is most strongly coupled is closed, the DgK* , but it also has sensitive couplings to 
all channels into which it is allowed to decay. The consideration of the finite width 
of the vector mesons has very small effect over this resonance (~2 MeV increase in 
its width). 

The couplings of the poles in this sector to the channels are given in Table [71 
As in the former cases, the states are clearly separated into the heavy vector-light 
pseudo-scalar and light vector-heavy pseudo-scalar sectors. 

3.4 C=1,S=-1,I=0 

The two remaining exotic members of the sextet should be found in this sector. One 
of them becomes a broad cusp at the KD* threshold when x = 1 while the other one 
is a narrow resonance with pole position (2756.08-12.15) MeV. The couplings of this 
pole are given in Tabled 
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Channel 


(2756.0^ 


^-i2.15) MeV 




\9i 


(GeV) 


DK* 


5.66 


KD* 


1.42 



Table 8: Residues for the C=1,S=-1,I=0 sector 



Channel 


(1101. 72-156.27) MeV 


(1213.20-10.89) MeV 




1^7.1 (GeV) 


1^7.1 (GeV) 


ttK* 


4.48 


0.51 


Kp 


1.57 


5.15 


r]K* 


0.36 


3.55 


Kuj 


3.02 


1.42 


DDI 


0.78 


0.16 


DsD* 


0.05 


0.48 


KJ/iP 


0.08 


0.02 


VcK* 


0.03 


0.02 



Table 9: Residues for the C=0,S=1,I=| sector 



When taking into account the 50 MeV width of the K* meson, this resonance gets 
a much bigger width, of the order of 180 MeV. In this case, and in all other sectors 
where the effect of the finite width of the vector mesons were taken into account, the 
only significant effect one could observe was over the width of the resonance. The 
mass of the resonances were affected in less than 0.5 % and the important couplings 
in less than 5%. 

3.5 C=0,S=1,I=| 

Two poles are found here coming from the two octets in the scattering of the low 
lying pseudo-scalar with the light vector mesons. In principle one could be tempted 
to assign these two poles to the two axial kaons from PDG [l3], but the mass of one 
of these, the iri(1400) is about 200-300 MeV off the pole positions we found and its 
width is much smaller than that. With this in mind we followed the interpretation 
of Roca |36j that the i^i(1270) should have a two pole structure. 

The couplings of the two poles to the different channels are in Table [9l 
This sector is explained in more detail in [36]. The novelty here is that, in spite 
of including now the heavy channels, the results are basically unaltered compared to 
those of [36] where only the light sector was used. This indicates a very weak mixing 
of the heavy and light sectors. 

Concerning the two Ki states it is also opportune to mention that in [16] some 
experimental information was reanalyzed giving strong support to the existence of 
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Channel 


(1230.15-i47.02) MeV 

h\ (GeV) 


(1012.95-i89.77) MeV 

h\ (GeV) 


vrp 




4.49 


KK* ± c.c. 


6.56 


2.42 


TTUJ 


3.07 




VP 


2.90 




DD* ± c.c. 


0.45 


0.88 


VcP 


0.03 






0.04 





Table 10: Residues for the C=0,S=0,I=1 sector 



these two states. 

3.6 C=0,S=0,I=1 

In this sector too there are two poles coming from the two octets but, since this is 
the non-strange sector, this two states have defined G-parity and therefore cannot 
mix. 

The pole with positive G-parity we associate with the 6i(1235) resonance. The 
small discrepancy between the experimental width and the value found from our 
theoretical model is explained since, experimentally, some decay channels of this 
resonance are three or four body decays while our model contemplates just two body 
hadronic decays. 

The negative G-parity pole should be identified with the ai(1260) but here the 
model gives a worse description of the resonance, the mass of the pole is smaller than 
expected although the huge width of the resonance makes this a minor problem. 
Also the width found within the model is very large, of the order of magnitude of 
the experimental one which is estimated with large errors. Again one should note 
that an important fraction of the width of this resonance could be due to many body 
decays not included in the present model. 

The couplings of the resonances to the channels are given in Table [TU] and they are 
very similar to those found in [36]. There the wg was substituted in terms of uj and 
(p and the sum of |gp for iruj and vr^ is similar to the \g\'^ for vrwg of our calculation. 

3.7 C=0,S=0,I=0 

Five poles are found in this sector. Tree have negative charge conjugation parity 
and two of them positive C-parity. In the light sector the positive C-parity pole is 
associated with the /i(1285), it appears in our model as a truly bound state, as it 
should, since none of its observed decay channels is a pseudo-scalar vector meson 
one, the possible decay channels within the model. The results obtained here and in 
the other two sectors for the light axial resonances without to — (j) mixing are very 
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similar to those obtained in [36] where the mixing was explicitly taken into account 
thus, corroborating the moderate effects of the mixing found for the charmed sector. 

The heavy singlet with positive C-parity obtained at 3837 MeV is a good candi- 
date to be associated with the controversial state X(3872). In this case this state 
is interpreted as being mainly a mixed molecule of DD* + c.c. and D^D* — c.c, its 
only possible decay channel within the model being the KK* + c.c. which is highly 
suppressed. In Table [H] the couplings of the two poles are presented. We can see 
there the strong decoupling of the heavy and light sectors. 

The low lying negative C-parity resonances can be associated with the two hi 
resonances. The singlet at (925. 12-124.61) MeV we identify with the /ii(1170) and, 
since we get it with a lower mass, our width is much smaller than the experimental 
one because our state has less phase space for decay. With the octet pole at (1213.00- 
15.67) MeV the same thing happens, and we associate it with the /ii(1380) despite 
the smaller mass and width compared with experimental values. 

In the heavy sector we find another state at 3840 MeV and negative C-parity. 
The decays of the X(3872), reported in [U], into ^J/ip and wJ/V' indicate that the 
C-parity of this state is positive. The decay into tt^tt^ J /ip reported in [2l] assuming 
that the ■k^tt~ comes from a p state would give the same C-parity but would imply 
isospin breaking. The large branching fraction 



S(X ^ vr+vr-J/V') 

reported in [U] indicates a massive violation of G-parity and hence isospin if one has 
only one X particle. 

There is a more appealing explanation for eq. (|27p if one had two X(3872) 
states with different G-parity and correspondingly C-parity. Should the tt+tt" in the 
denominator of eq. ()27p correspond to an 1=0 state one would not have to invoke 
isospin violation, but instead the existence of a negative G-parity (and hence C- 
parity) state. This would imply that there is strength of these events in the a region 
of the TTTT invariant mass, and this seems to be the case as reported in [27], although 
the statistics is low. This latter scenario would fit with our predictions of two states 
nearly degenerate with opposite C-parity. A similar argumentation has been made in 
|48j to justify the existence of two degenerate X(3872) states with different C-parity. 

The difference of 35 MeV in the binding energy between our model and experiment 
(~ 1% difference) is perfectly acceptable for a theoretical model that looks at the 
whole spectrum of axial vector mesons in a broad range of masses with only two 
free parameters {an, o-l)- Yet, we can do some fine tuning to get a mass like in 
experiment by changing the subtraction constants in the loops. In this case we can 
take an = —1.30 (from —1.55 before) and then we find the mass of the positive 
C-parity pole at 3872.67 MeV and simultaneously the state with negative C-parity 
just disappears as a pole, by crossing the DD* threshold, and leads to a marked 
cusp structure in the DD* — c.c. amplitude. Indeed, changing an to slightly more 
negative values we regain the pole just below this threshold. 

Recently a new possible state at 3875 MeV was reported at Belle [50] decaying into 
DqDq. The 3 MeV difference of this state with the X(3872) is precisely the difference 
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E{MeV) E(MeV) 



Figure 2: a)|Tp for the positive G-parity state in the DD* channel compared with the 
Belle Data (in this plot uh = —1.23) b) |Tpp for both G-parity states in the DD* channel 
compared with the Belle Data (in this plot au = —1.30 for the G=- state) 



of masses between the positive and negative C-parity states that we obtain. It is thus 
tempting to associate to the new X(3875) state our negative C-parity state. Such 
scenario is not excluded by the data as we show in the analysis below. 

In the decay of a B particle to K7r~^7r~ J/ip we have, defining E = Minv{TT^ J /ip) — 



^^''+ oc |T(£»OZ)°* + c.c. ^D°D0*+c.c.)pr(X^7r+^-J/V') (28) 



clE 



which is approximately proportional to |Tp in a few MeV range of E around M£)0£)0* 
since T{X — > 7r~^7r~J/ip) barely changes in this range of energies given the relatively 
large phase space for this decay. On the other hand for the B particle decaying into 
KD^D^* we have (now E = Mi„„(Z)°L>°*) - M^odo*): 



oc \T{D^D^* - c.c. ^ D^D^* - c.c.)\'^T{X' ^ D^DO*) (29) 

uE 

which is approximately proportional to where p is the three-momentum of the 

D particle in the center of mass frame of the DD* system, since now T{X' D^D^*) 
is proportional to this three-momentum. 

In fig [2] we plot -^g^ and -^g- as a function of E and we compare our results with 
the experimental data from Belle [50] [51]. We can see that in the case of the positive 
C-parity state the tt^tt" J/-^ distribution is very sharp, while the distribution of the 
invariant mass of D^D^* leads to a strong enhancement of the distribution around 
the D^D^* threshold, in both cases in fair agreement with experiment. 

We should note that an empirical analysis of the data in a recent paper [52] , taking 
only one resonance, produces a similar behavior assuming that the resonance couples 
very strongly to D^D^*, as it is also our case (see couplings in Tables [TTHl2p . In view 
of the results of [52] it is also instructive to see what our model would give assuming 
that the resonance with positive C-parity is responsible for the two distributions in 
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Channel 


1292.96 MeV 


(3837.57-iO.OO) MeV 




\g,\ (GeV) 


\g,\ (GeV) 


^^(DD* + c.c.) 


0.15 


13.61 


^^{D,D: - c.c.) 


0.54 


10.58 


^^{KK* + c.c.) 


7.15 


0.03 



Table 11: Residues for the C=0, 8=0,1^=0"'" sector 



Channel 


(925. 12-124.61) MeV 


(1213.00-15.67) MeV 


(3840.69-11.60) MeV 




\g,\ (GeV) 


\g,\ (GeV) 


\9:\ (GeV) 


Tip 


3.65 


1.07 


0.03 




0.10 


4.19 


0.00 


:kiDD*-c.c.) 


3.93 


1.03 


13.44 




0.00 


0.04 


1.25 


TjJ /ip 


0.01 


0.06 


1.20 


-1={DsD: + c.c.) 


2.31 


1.26 


9.96 


^,{KK*-c.c.) 


0.73 


6.14 


0.01 




1.08 


0.34 


1.99 



Table 12: Residues for the C=0,S=0,F=0- sector 



fig. [21 We show the result obtained for \T\'^p with the positive C-parity resonance in 
figure E] (b) with dotted line. As we can see, the shape and strength obtained in this 
case is not very different from that of the negative C-parity resonance. Hence, this 
alternative scenario, which would correspond to the one in [52], is not ruled out by 
these combined data. 

As mention above the ratio of eq. ()27p and the possible strength seen for the 
B{X 'K^'K~ J/ip) in the region of invariant vrvr masses around rria-, is so far the 
strongest indications favoring the two C-parity states. 

4 Overview and Conclusion 

We studied the dynamical generation of axial resonances by looking for the poles in 
the scattering T-matrix of pseudo-scalars with vector mesons. For the interaction 
Lagrangian we first constructed a 5^7(4) flavor symmetrical Lagrangian for the in- 
teraction of the 15-plet of pseudo-scalar mesons with the 15-plet of vector mesons. 
The symmetry was broken down to SU{'i) by suppressing exchanges of heavy vector 
mesons in the implicit Weinberg- Tomozawa term, following a prescription developed 
in a previous paper [34]. From the Lagrangian, tree level amplitudes were evaluated, 
projected in s-wave and collected in a matrix for the various channels. This matrix 
was used as the potential, transformed to an isospin basis, in the Bethe-Salpeter 
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equation, which provides the unitarized ampUtudes between the channels. 

The poles generated within the model can be associated with the various axial 
resonances listed by the Particle Data Group [l3] , and also many new resonances are 
predicted: three broad ones in the mass range between 2.5 and 2.6 GeV, and three 
narrower ones around 2.75 GeV, these resonances belong to two SU{3) sextets, they 
have, therefore, exotic quantum numbers and they have not yet been experimentally 
observed. We summarize these states with their quantum numbers in Table [T3l 



c 


s 


IG^JPC^ 


Re{^) (MeV) 


lra{^) (MeV) 


Channel 


1 


1 


1(1+) 


2529.30 


-238.56 


nD*, KD* 


1 







Cusp (2607) 


Broad 




1 


-1 


0(1+) 


Cusp (2503) 


Broad 




1 


1 


1(1+) 


2756.52 


-32.95 [cusp] 


DsP, DK* 


1 





1(1+) 


2750.22 


-99.91 [-101] 


Duj 


1 


-1 


0(1+) 


2756.08 


-2.15 [-92] 


DK* 








o-(i+-) 


3840.69 


-1.60 





Table 13: List of predicted and not yet observed resonances with quantum numbers and 
the open channels to which they couple most strongly. The results in brackets for the 
I■m^/s are obtained taking into account the finite width of the p and K* mesons. 



Remaining discrepancies between our model and experiment can be attributed 
to possible many body decays of these objects or possible higher order terms that 
could be included in the Lagrangian. The size of the discrepancies, typical of any 
successful hadronic model describing hadronic spectra, can be used to estimate the 
uncertainties in the predictions for the new states that we obtain. 

Some states obtained here have been also studied before in a similar framework, 
but with more restricted coupled channels space. The similarity of the results rein- 
forces these findings. Many other states are reported here for the first time within this 
unitary coupled channel framework. In the light scalar sector the poles and couplings 
found within our approach coincide very well with the ones found by Roca [36] . This 
happens because, despite the enlarged coupled channel space including heavy mesons, 
this new sector couples very weakly with the light one. In the open-charm sector the 
poles found for the lightest anti-triplet coincide with the results found by Kolomeit- 
sev [9] and Guo [TT]. As already happened for the scalar resonances, the poles found 
within our model for the sextet state in this sector have broader widths because of 
the use, in our model, of a different meson decay constant for heavy mesons. Besides, 
our model allows also for the inclusion of channels with heavy pseudo-scalar mesons 
interacting with light vector ones, as a result of which our model generates a richer 
spectrum, with poles for an extra anti-triplet and an extra sextet. In the charmed 
sector, C=l, we find six resonances not yet observed. Three of them are either too 
broad, or they degenerate into cusps as SU{3) is gradually broken. However, three 
of them remain sufficiently narrow, such that they could in principle be detected. 

Moreover our Lagrangian incorporates the hidden-charm sector and an attractive 
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interaction in the 3 (Xi 3* and 3 ® 3* charmed mesons is responsible for the generation 
of two resonances. One of them can be associated with the new X(3872) state. The 
other one appears as a strong cusp in our case, or a resonance with shghtly increased 
attraction, and could be associated to the peak structure found recently around the 
D^D^* threshold at Belle, although this structure can also be explained in terms of 
only one state. On the other hand, the ratio of partial decay width into two or three 
pions and J ftp (eq. (I27p ) provides strong support for the existence of two states with 
different C-parity. 

The agreement of the theory with data for the known resonances, together with 
the success of the theory for the charmed scalar mesons [3l] , gives us confidence on 
these new predicted states, such as to strongly suggest their experimental search. 
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A The ^ Coefficients 



A.l Open-Charm Sector (C=l) 



s= 


=1,1^ 




:1(1+) 




Channels 


^Dl 


DsP 


KD* 


DK* 


ttD* 








-1 


-7 


DsP 








-7 


-1 


KD* 


-1 


-7 








DK* 


-7 


-1 









S=L I'''(./'''")=0(1+) 



Channels 



DK* 



KD* 



riDt 



DsU 



DK* 
KD* 
vD*s 

DgLO 

VcD*s 
DJ/^ 



-2 
-2 

VS 



7 

-V3 ^ 



-2 /2 



V3 

§7 



7 

V3 



27 
3 



2^17 



V3 -2 ^§7 




7 

V3 
27 
3 



2V27 
3 



2V27 
3 



47 
3 



2V27 
3 



47 
3 



S=0, I«(J^^)=i(l- 



Channels 


ttD* 




KDt 




-2 


7 
2 




Dp 


7 
2 


-2 





KDl 


fs 
V 2 





-1 


DsK* 





/? 

V 2 





r]D 





7 
2 


V 2 


Du 


7 
2 





-Vl7 


ricD* 





V27 





DJ/^P 


-V27 





2l 
V3 



DgET* r/Z) Dlv r]cD* DJ/ip 







§7 



27 






7 
2 



3 
2 

§7 





7 



3 



§7 



727 
3 





\/27 



27 
Vs 


V27 

3 



47 
3 



-\/27 



27 

%/3 



V27 

3 



47 
3 



S=-l, I«( 


J^^) = 


=0(1+) 


Channels 


DK* 


KD* 




-1 


-7 


KD* 


-7 


- 1 
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Cliamiels 



A.2 Hidden-Charm Sector (C=0) 

S=l, I«(J^^)=i(l+) 

Kp 



riK 



DDI D,D* KJ/y ijcK" 



ttK* 
Kp 

rjK* 
Ku 

DDI 

DsD* 
KJ/il; 

VcK* 



-2 


1 

2 





3 
2 


1 

2 


-2 


3 
"2 








3 
2 





3 
2 


3 
2 





3 
2 





\/i7 





7 


^7 







717 


7 

V6 



























§7 







7 

§7 





27 
V3 







h 

h 



7 



27 
V3 



























27 


27 


V3 


V3 


27 


27 


V3 


V3 















S=0, lg(J^^)=l+(l+-) 



Channels 


TTiO 


VP 


-^^{KK* + c.c.) 


-^{DD*+c.c.) 


VcP 


TrJ/tp 


TTiO 








V3 


7 

V3 








VP 








V3 


7 

V3 








j.{KK* + c.c.) 


V3 


V3 


-1 


-7 








j=,{DD* + c.c.) 


7 

73 


1 

V3 


-7 


-V 






VcP 











2^17 




























S=0, lg(J^g)=l-(l++) 

Channels | - c.c.) 7rp -j^{DD* - c.c.) 

^{KK*-c.c.) -1 -V2 7 

Trp -V2 -2 a/27 

^(5L»*-c.c.) 7 V27 -V' 



S=0, I^(J^^)=0+(1++) 

Channels \^{KK* + c.c.) ^{DD*+c.c.) ^{DsD* - c.c.) 
-^{KK* + c.c.) Is ^ 72^ 

hOsD^-cc.) V27 -V2 -V-l 
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s=o, ig(j^g)=o-(i+-) 



Channels 


np 


rjuj 


—c.c.) 




r]J/ip 


+ C.C.) 


—c.c.) 


rjcJ/ip 


" P 


-4 





a/37 











V 













7 
3 








2v^7 
3 


3 





j^^{DD*-c.c.) 




7 
3 


-ip-2 


2V27 
3 


2V27 

3 


-V2 


7 


87 

3 


rjcUJ 








2V27 
3 








47 
3 
















2v^7 
3 








47 
3 








-^{DsDt+cc.) 





2^27 
3 


-V2 


47 
3 


47 
3 


-tp-l 




4^27 

3 


j=^{KK*-c.c.) 


-v^ 


3 


7 








\/27 


-3 













87 
3 








4V27 
3 
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